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^ SUMMARY 

Flight  data  from  a 60f,  delta  wing  aircraft  have  been  analysed  using  a modified 
Newton- Raphson  parameter  estimation  procedure.  The  model  equations  used  for  the 
analysis  were  extended  to  account  for  incidence  vane  errors  and  non-linearities  in  the 
pitching  moment  curves.  Longitudinal  derivatives  extracted  from  the  data  have  been 
compared  with  wind  tunnel  measurements  and  some  theoretical  estimates  and  areas  of 
agreement  and  disagreement  identified.  The  results  demonstrate  the  usefulness  of  the 
parameter  identification  method  not  only  for  the  validation  of  aircraft  mathematical 
models  and  for  checking  flight  results  against  wind  tunnel  data  but  also  for  obtaining 
aerodynamic  data  not  easily  available  through  other  means. 
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I.  INTRODUCTION 


The  use  of  systems  identification  techniques  has  been  developing  into  a standard  method 
for  flight  test  analysis  over  the  past  few  years.  They  provide  for  reduced  flight  testing  time  and 
more  comprehensive  data  analysis*  and  have  been  applied  to  a variety  of  problems  including 
high  angle  of  attack  “wing-rock”  analysis^  and  determination  of  in-flight  airload  parameters® 
as  well  as  the  more  common  types  of  flight  manoeuvre.  From  the  point  of  view  of  mathematical 
modelling  they  provide  a means  for  rational  and  systematic  validation  of  a model  and  a means 
for  obtaining  additional  data,  which  may  not  be  available  from  wind  tunnel  tests,  for  use  in 
such  models. 

A summary  of  various  parameter  estimation  procedures  applied  to  aircraft  can  be  found  in 
Reference  4,  for  example.  One  of  these  methods,  the  Modified  Newton- Raphson  technique  as 
described  in  References  5 and  6,  has  been  adapted  for  use  on  the  ARL  PDPIO  computer.  It  was 
studied  in  some  detail  in  Reference  7 with  a view  to  using  it  in  support  of  mathematical  modelling 
activities.  Briefly,  referring  to  Figure  I,  the  method  sets  out  to  minimise  a cost  functional  J, 
proportional  to  the  difference  between  the  measured  response,  z,  and  the  calculated  response,  y, 

J = (I/A')  2 (Zi  - yO^^>i(*<  - y<)  + (c  - CoVDi{c  - Co)  (I) 

!=  1 

where  Di  and  Dz  are  diagonal  weighting  matrices  and  the  second  term  is  a weighted  mean  square 
difference  between  the  parameter  vector,  c,  and  a stipulated  a priori  value,  Co. 

The  mathematical  model  is  described  by  the  state  equation 


iL^  A\  f Bu 


and  the  calculated  response  is 

y = Fx  + Cu  -{  • b 

The  measured  response  is  assumed  to  be  of  the  form 


where 


z = y -f-  n 


(2) 

0) 

(4) 


X is  the  vector  of  state  variables 
u is  the  vector  of  control  variables 
y is  the  response  vector 
b is  a measurement  bias  vector 
z is  the  measured  response  vector 
n is  a random  noise  vector 

A,  B,  Fand  G are  matrices  of  parameters  defining  the  model. 

The  parameter  vector,  c,  contains  some  or  all  of  the  elements  of  A,  B,  Fand  G,  the  elements 
of  b and  the  initial  conditions.  The  iterative  Modified  Newton-Raphson  procedure  is  used  to 
find  the  value  of  c which  minimizes  the  cost,  J.  An  indication  of  the  quality  of  the  estimates,  c, 
is  given  by  the  Cramer-Rao  bound,  ocr,  which  is  a lower  bound  on  the  variance  of  the  estimates.** 
While  Reference  7 considered  more  general  aspects  of  the  parameter  estimation  process 
including  questions  of  response  information  requirements,  control  input  design,  data  sampling 
rates,  etc.,  this  note  describes  the  application  of  the  Modified  Newton-Raphson  technique  to 
the  analysis  of  actual  flight  data  in  the  form  of  aircraft  longitudinal  response  to  elevator  input. 
The  detailed  model  equations  are  outlined  in  Section  2.  These  include  allowance  for  incidence 
vane  errors  on  the  one  hand  and  extension  of  the  model  to  treat  aerodynamic  non-linearities  on 
the  other.  Section  3 describes  the  flight  data  including  data  conditioning  prior  to  use  in  the 
analysis,  and  discusses  the  a priori  values  and  weights  used  in  the  matching  process.  The  results 
of  the  analysis  and  general  discussion  follow  in  Section  4 which  also  includes  some  comparison 
of  the  values  of  the  extracted  parameters  with  wind  tjnnel  and/or  theoretical  values. 
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2.  MATHEMATICAL  MODEL 

In  this  section  the  detailed  form  of  the  mathematical  model  described  in  general  by 
Equations  2 and  3,  is  developed.  Firstly,  the  basic  linear  small  perturbation  equations  for  aircraft 
longitudinal  response  are  outlined.  Next  a simple  extension  of  the  equations  is  made  in  order 
to  treat  aerodynamic  force  and/or  moment  non-linearities  within  the  context  of  a linear  model. 
A further  extension  of  the  model  then  follows  in  order  to  deal  with  errors  in  the  incidence  vane 
measurements.  Finally,  one  method  for  accommodating  instrument  lags  and/or  calibration 
errors  is  described. 

2.1  Basic  Model 

It  is  assumed  that  the  aircraft  is  trimmed  in  steady  level  flight  with  velocity  K,,,  incidence 
and  attitude  Of.  Small  perturbation  motions  about  this  initial  state  are  produced  by  elevator 
inputs.  The  linearised,  uncoupled,  longitudinal  equations,  written  in  body  axes,  to  describe  this 


motion,  are  (see  for  example  Reference  8 for  more  details); 

mu  = XuU  + -f  X»ie  — mw,q  — mg  cos  (5a) 

mw  = ZuU  f Z^a  I (Z»  f mue)q  '•  ZJi  f Z,h,  — mg  sin  M (5h) 

lyq  = A/uM  f f Mji  f Mfht  (5c) 

6 = q (5d) 

where  We  = sin  ««  ~ and  Ue  = F,  cos  ««  ~ K,  for  small 


For  the  short  period  mode  it  is  further  assumed  that  airspeed  remains  constant  at  the  trim 
value,  V = Ve-  Then,  from  the  definitions 

(w  f W,)  = Ve  sin  (a  4 a,)  and  (u  f U,)  = Ve  cos  (a  + Xe) 
it  follows  that  (for  small  x and  ««); 

w = VeX  (6) 

and  u = —Ve(x,  0 -5x).  (7) 

Thus  Equation  5b  for  w becomes  an  equation  for  x via  Equation  6 while  Equation  5a  for  ii 
becomes  redundant  since  u is  related  to  « through  Equation  7.  Substituting  for  u from  Equation  7 
into  Equations  5b  and  5c  leads  to  the  short  period  approximations; 

(mV,  — Zi)dc  = [Z,  — Ke(o,  4-  0-5a)  Zu]x  f (Z,  4 mue)q  t Ztb,  — mg  sin  6,0 


~ Z^x  4-  mueq  + Z,b,  — mg  sin  M (8fl) 

lyq  = [A/a  — Ve(xe  f 0-5a)A/«]a  4-  A/^  4“  M jjx  4 M 

~ A/a«  f A/^  4 A/ii  4 M,he  (86) 

6 =^q.  (8c) 

These  equations  can  be  further  simplified  by  noting  that  Z^  mV,.  Substituting  for  x in 
Equation  86  from  Equation  8a  gives 

lyq  = (Ma  + MiZJmVe)x  + (A/,  4-  M^UelVe)q  4-  (A/,,  4-  MiZtJmVe)^,  - Mig  sin  M/K, 

=;  M,x  + (My  f Mi)q  4-  M,he  (8rf) 

where  it  has  been  assumed  that  u,  ~ V,  and  the  neglected  terms  are  either  much  smaller  than 


those  retained  or  else  make  very  little  difference  in  the  parameter  estimation  calculations.  In  each 
case  the  assumptions  have  been  justified  by  the  subsequent  results. 

An  equation  for  the  normal  acceleration,  nt,  follows  from  the  definition 

gmlVe  = (w  — u,q  + g sin  6,6)1  V,  ~ (ZJmV,)x  4-  (ZylmVt)q  4-  (Z,JmVe)Se-  (9) 

Written  in  matrix  notation  the  mathematical  model  used  to  analyse  the  short  period  motion 
can  be  summarized  as  follows.  The  state  equations  become; 
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it 

'ZJmVe  1 

—g  sin  0,1  Vi 

a 

'ZtJmV, 

z„ 

4 

= 

MJIy  (My  1 Mi)/Iy 

0 

Q 

1 

M,Jly 

My 

.1. 

6_ 

. 0 1 

0 

0 

0 

0 

The  equations  for  the  response  variables  become ; 


a 

" 1 0 0 

a 

0 

0 

S* 

Q 

0 1 0 

9 

1 

0 

0 

1 

0 

0 0 1 

0 

■r 

0 

0 

gntlVe, 

ZJmV,  Zy/mV,  0_ 

ZiJmV,  Zo_ 

The  parameter  vector  is  given  by: 

c = [ZJmVe,  ZsJmVe,  MJIy,  (My  | Mi)/Iy,  M,JIy,  ZylmVe,  Zo, 


(10) 


(11) 


(12) 


where  Zo  and  Mg  have  been  introduced  to  account  for  possible  bias  in  the  state  and  the 
measurement  bias,  b,  appearing  in  Equation  11,  has  been  left  out  of  the  parameter  vector  as  it 
will  be  neglected  from  here  on. 


2.2  Model  with  Non-linear  Aerodynamics 

The  basic  model  of  Section  2.1  assumes  linear  aerodynamics.  However  it  was  obvious  from 
inspection  of  the  wind  tunnel  results  of  Reference  9 that  at  transonic  speeds  the  pitching  moment 
curve  can  be  expected  to  be  significantly  non-linear  with  incidence,  a.  Over  the  present  range  of 
interest,  i.e.  2°  < a a«  < 8°,  the  wind  tunnel  curves  suggest  that  the  moment  increment  due 
to  incidence  needs  to  be  modelled  by  a cubic  equation,  i.e. 

AM  = M^a  1 M^ia^  + MjjaS.  (13) 

In  order  to  accommodate  this  within  the  framework  of  a linear  model,  the  ofl  and  a®  terms 
can  be  treated  as  known  inputs  with  a®  and  a®  becoming  part  of  the  input  vector  u together 
with  8f.  The  time  histories  of  ofl  and  a®  may  be  calculated  directly  from  the  measured  values  of 
incidence,  a^,  or  alternatively  a purely  linear  model  may  be  used  initially  with  the  resulting  a 
time  histories  used  to  provide  a®  and  a®  for  a subsequent  matching  in  an  iterative  manner. 
As  long  as  the  non-linearities  are  relatively  small  this  approach  appears  to  be  successful  and  the 
two  alternatives  for  obtaining  the  h ories  of  a®  and  a®  produce  only  small  differences  in  the 
final  results. 

The  only  equation  altered  by  the  inclusion  of  the  “non-linear”  terms  is  Equation  8</  for 
pitch  acceleration  q,  which  now  becomes; 


lyq  = Mg<x  1 (M,  -I-  Mi)q  I M,^S,  + M^^a®  i M„,a®.  (14) 

The  matrix  equations  (Equations  10  and  II)  are  modified  slightly  with  the  control  vector 
now  becoming; 

[Se,a2,a3,  (15) 

and  matrices  B,  G become  4x4  and  5x4  matrices  respectively  with  the  additional  elements 
being  zero  with  the  exception  of  the  two  terms  in  B which  correspond  to  M^j/Zy  and  M^i/Iy. 
Finally  the  parameter  vector,  c (Equation  12),  is  increased  by  the  inclusion  of  M^illy  and  M^yjly 
which  need  to  be  identified. 

In  some  cases  the  Z-force  also  appears  to  be  a non-linear  function  of  incidence  and  a similar 
approach  to  the  above  may  also  be  applied  there.  However,  unless  it  is  thought  important  to  the 
model  to  include  such  additional  parameters,  the  identification  process  may  often  be  carried 
through  quite  successfully  without  them.  That  is  to  say,  while  the  addition  of  such  terms  may 
make  an  improvement  to  the  time  history  matches  achieved,  the  estimated  values  of  the  other 
parameters  of  interest  may  differ  very  little.  At  the  same  time  computational  effort  would  increase. 
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2.3  Corrections  to  Incidence  Vane  Measurements 

The  incidence  measurements  are  influenced  by  a pilch  rate  contribution,  a,,  described  as 
follows: 

«/»  = « — «,  (16) 
where  «,  = (xviyt)q 

and  Xv  is  the  distance  of  the  vane  ahead  of  the  centre  of  gravity. 

Equation  16  can  be  modified  to  account  for  possible  vane  calibration  slope  errors  as  follows: 

*m  = G a — a,  (17) 

where  C has  to  be  estimated  from  the  matching. 

The  vane  was  not  dynamically  balanced  so  that  the  strain  gauges  also  responded  to  vane 
accelerations  (see  Reference  10  for  more  details).  The  correction  to  measured  incidence  due  to 
normal  acceleration  at  the  vane  can  be  written  : 


ad  = K lq(x„/  ye)  - gnil  ye]  (18) 

where  pitch  acceleration,  q is  in  radians/sec^  and  normal  acceleration,  is  in  g’s,  while  the 
constant  K depends  on  vane  characteristics.  The  measured  incidence,  am,  can  now  be  related  to 
the  true  incidence,  a,  by  the  equation : 

am  = Ca  — a,  — otd  (19) 

= Ga  - (Xv/ye)q  - Klq(Xe/ye)  - gHe/ye] 

or  am  + *8  = Ga  — K[q{Xviye)  — (20) 

By  substituting  for  and  q from  Equations  9 and  14,  Equation  20  can  be  expanded  to  give: 

«m  + a,  = {G  - K[(Xv/ye)  MJIy  - Zjmye]}  a - (A/,  -p  - 

- K[(x^lye)  M,Jly  - Z,Jmye]&e  - K(Xeiye)  {M,,IIy)ofl  - (Afa^//,,)*^ 

= flia  -t-  029  + bihe  + biofi  -f  63a®.  (21) 

The  correction  due  to  pitch  rate,  a,,  can  be  immediately  calculated  from  Equation  16  since 
measurements  of  q are  available.  Thus  (am  -V  a,)  rather  than  am  is  treated  as  a response 
variable  which  can  be  compared  with  the  model  results  as  calculated  from  the  right  hand  side 
of  Equation  21. 

The  full  system  of  equations  including  non-linear  pitching  moment  model  and  incidence 
vane  corrections  is  summarized  below. 

The  state  equations  are: 


a 

q 

— 

J 

Zjmye  1 -gsin0e/ye 

MJIy  (My  + MJlIy  0 
0 1 0 


The  response  equations  are: 

I 0 0 

0 1 0 

0 0 I 

01  a-i  0 


a 

q 

e 

= 

am  + a, 

- g^ijye . 

a 

'Z,Jmye 

0 

q 

■f 

M,JIy 

MoJIy 

0 

0 

a 

■ 0 

0 0 

q 

} 

0 

0 0 

1 9 

0 

0 0 

bi 

bt  bt 

Z,Jmye 

0 0 

Se- 

a® 

. 1 _ 

0 

0 

0 

bo 

Zo 


Izjmy,  ZJmy.  0 
The  parameter  vector  now  is : 

c = [Zjmye,  Z,Jmye,  MJIy,  (My  MJ/Iy,  M,Jly,  MJ/y. 

Meiijly,  zjmye,  a\,  at,  b\,  bt,  bt,  Z,,  Mo, 


aZ 

as 

_ 1 . 

(22) 


(23) 


(24) 
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where  bo  has  been  added  to  account  for  possible  state  bias.  The  measurement  bias  vector,  b, 
has  again  been  assumed  to  be  zero.  The  coefficients  ai,  Oi,  bu  bz  and  bz  have  been  treated  as 
unknown  parameters  to  be  estimated.  Alternatively,  an  estimate  for  K can  be  made  from  the 
known  inertial  and  aerodynamic  characteristics  of  the  vane  and  approximate  values  of  these 
coefficients  established.  The  less  critical  ones,  such  as  bz  and  bz  may  be  fixed  and  only  ai,  az 
and  bi  then  allowed  to  vary.  The  subsequently  identified  values  of  az  and  b\  can  then  be  used 
in  the  definitions  of  az  and  b\  (Equation  21)  to  obtain  two  values  for  K.  Hopefully  these  values 
will  be  in  reasonable  agreement  with  each  other  and  with  the  previously  estimated  value. 
Having  obtained  K in  this  way  a value  for  G follows  from  ai  (Equation  21). 


2.4  Instrument  Lags  and  Calibration  Errors 

Calibration  errors  and  lags  in  instrumentation  can  be  modelled  in  the  state  equations. 
For  example,  if  a lag  and/or  calibration  change  is  suspected  in  the  incidence  measurements, 
the  following  model  could  be  postulated : 


, am  G(x 

dt  T T 


(25) 


Where  r is  the  lag  time  constant  and  G would  account  for  any  change  in  calibration  slope. 
The  values  of  t and  G would  then  be  identified  for  optimum  match.  Such  an  approach  was 
initially  used  to  improve  the  match  of  the  incidence  record  and  although  an  improvement  was 
achieved,  as  may  be  expected  when  two  additional  free  parameters  are  provided,  it  was  obvious 
from  the  results  that  this  approach  was  inadequate  and  much  better  results  were  obtained  using 
the  corrections  outlined  in  Equation  2.3.  This  will  be  more  fully  discussed  in  Section  4. 

3.  FLIGHT  DATA 

This  section  discusses  the  flight  data  to  be  analysed.  The  form  of  the  data  and  the  conditions 
under  which  they  were  obtained  are  first  outlined  followed  by  a description  of  corrections  and 
pre-processing  applied  to  the  data  in  preparation  for  the  Modified  Newton-Raphson  program. 
Finally,  the  a priori  values  used  for  each  of  the  parameters  are  summarized  followed  by  a 
discussion  of  the  values  chosen  for  the  elements  of  the  weighting  matrices  Di  and  Dz- 

.3.1  Test  Conditions 

The  data  available  relate  to  flight  tests  performed  on  a 60°  delta  wing  fighter  aircraft. 
The  manoeuvres  considered  in  the  present  study  are  aircraft  pitch  response  to  pilot  input  in  the 
form  of  an  elevator  pulse.  Flights  at  two  different  Mach  numbers  were  analysed,  namely 
M = 0-96  and  0-71.  Other  relevant  conditions  for  these  flights,  referred  to  as  flights  210  and 
211  respectively,  are  summarized  in  Table  I. 


TABLE  1 

Flight  Test  Conditions 


Details 

I 

Flight  210 

Flight  211 

M 

096 

0-71 

h 

9019  m (29  590  ft) 

2971  m (9746  ft) 

Ve 

292 -2  m/s  (568  kt) 

233  -8  m/s  (454  kt) 

«« 

3-64° 

2-64° 

f>e 

3-28° 

2-47° 

(S«). 

-0-63° 

-0-82° 

eg- 

50  -49%  of  chord 

49-10%  of  chord 

K 

0-293 

0-234 

The  values  of  K given  in  Table  1 were  estimated  from  data  given  in  Reference  10.  From 
Equation  18  it  follows  that 

KgjV,  = (idl[q{Xvlg)  - ni],  rad/g 

Now  for  an  effective  vane  mass  of  0-325  kg,  the  inertial  force  per  g = 318  N.  On  the 
other  hand,  from  wind  tunnel  calibrations,*®  force  per  radian  = 323-074  N so  that  equivalent 
radians  per  g for  the  inertial  force  = 3- 18/323-074  = 0-00984  = KgjVe.  Thus 

/f  = 0-00984  (26) 

Measured  data  consisted  of  response  time  histories  of  a,  q,  rii  and  8^.  Approximately  15 
seconds  of  record  was  available  but  only  the  5 seconds  following  (and  including)  the  pulse  input 
was  used  in  the  analysis.  This  represented  approximately  2 to  3 cycles  of  oscillation. 

The  data  rate  was  60  samples/second  and  from  inspection  of  the  flight  records  it  was 
determined  that  the  natural  frequency  w„  for  flight  210  was  approximately  3-9  rad/s  while  for 
flight  21 1 it  was  approximately  4-5  rad/s.  Using  an  in-between  value  of  4-2  rad/s  it  follows  that 

a>„A/  = 4-2/60  =0-07 

4-2  X 5 = 21  (27) 

where  At  is  the  sampling  interval  and  T is  the  total  record  length. 

The  values  of  to„A/  and  <unT  from  Equation  27  compare  with  the  values  of  conAr  < 0-14 
and  w„T  > 14  found  in  Reference  7 to  be  adequate  for  matching  of  the  Dutch  Roll  response 
to  rudder  doublet  input.  In  both  cases  the  systems  are  predominantly  second  order  though  the 
present  case  differs  from  that  of  Reference  7 in  having  a pulse  rather  than  a doublet  input. 
On  the  other  hand  the  present  system  of  equations  is  simpler  with  fewer  parameters.  It  may  be 
opected  then  that  the  sampling  rate  and  record  length  requirements  for  successful  identification 
'vould  be  similar  and  in  fact  no  difficulties  of  convergence  or  identification  were  encountered  in 
the  present  analysis.  Less  stringent  values  of  conAr  and  wnT  than  those  in  Equation  27  would 
also  be  expected  to  lead  to  successful  analysis.  Although  10  iterations  were  normally  specified 
in  the  program,  good  convergence  was  generally  achieved  by  about  5 iterations  or  less  with  only 
small  improvements  subsequently. 

3.2  Data  Conditioning 

Before  proceeding  with  the  analysis  a certain  amount  of  manipulation  and/or  correction  of 
•Jata  was  necessary.  This  included  the  following: 

(a)  Normal  accelerometer  corrections  for  location  errors  according  to  the  equation 
(Reference  8); 

Hz„  — (Xaq  — yap  Zaq'^  -\  ZaP^)lg  (28) 

where  is  the  measured  acceleration,  is  the  value  at  the  centre  of  gravity  and 
Xa,  ya,  Za  are  the  co-ordinates  of  the  accelerometer  relative  to  the  aircraft  body  axes. 
If  measurements  of  angular  accelerations,  p and  q,  are  not  available,  the  identification 
can  proceed  initially  without  this  correction  and  if  they  are  thought  to  be  important 
the  values  calculated  via  the  identified  model  used  in  a second  iteration.  In  the  present 
work,  these  corrections  were  generally  small. 

(M  Pitch  rate  gyro  corrections  for  misalignment  according  to  the  equations  given  in 
Reference  8 (p.  31).  Using  the  gyro  alignment  data  the  corrected  pitch  rate,  q,  is  obtained 
from  the  measured  roll,  pm,  pitch,  ^m,  and  yaw,  r^,  rates  according  to  the  equation 

? = 0-0302 />m  f I -00087  -0-0082  r™.  (29) 

(<')  Angle  of  incidence  corrections  for  pitch  rate  effect  as  discussed  in  Section  2.3.  The 
corrected  incidence  becomes  am  f ««  where  am  is  the  measured  value  and  a,  is  given 
by  Equation  16. 

(</)  A smooth  curve  was  fitted  by  hand  to  the  elevator  input  data  as  shown  in  Figures  2 and  3. 
The  effect  of  this  smoothing  on  the  results  was  hardly  discernible,  especially  for  flight  210 
where  the  measurement  noise  was  small  but  for  cases  with  higher  noise  levels  on  the 
measured  input  smoothing  should  minimize  the  standard  deviations  of  the  identified 
parameters  and  should  be  worthwhile  provider'  inadvertent  errors  are  not  introduced  in 
the  process. 
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(e)  Non-zero  steady  state  values  of  elevator  angle,  (S,),,  incidence,  and  normal  acceler- 
ation, (=1),  were  subtracted  from  the  respective  records  so  that  the  increments 

remaining  could  be  compared  directly  with  results  of  the  small  perturbation  model. 

3.3  A Priori  Values  and  Weighting  Matrices 

The  a priori  values  which  were  ascribed  to  the  various  parameters  were  mainly  obtained 
from  wind  tunnel  tests  (e.g.  Reference  9)  or  else  were  estimated  from  data  sheet  sources.  Table  2 
summarises  the  values  used  in  the  two  cases  analysed. 


TABLE  2 

A Priori  Parameter  Values 


Parameter 

Flight  210 

Flight  211 

ZJmVe 

-0-848 

-1-145 

ZsJmVe 

-0-314 

-0-467 

MJ/y 

-25-49 

-19-78 

(A/,  4 Mi)IIy 

-1-02 

-2-38 

M,JIy 

-35-46 

-31-14 

Maolly 

0 

— 

MJIy 

0 

— 

ZylmVe 

0 

0 

a\ 

0-93 

0-87 

aa 

0-0091 

0-0134 

bi 

0-24 

0-231 

bz 

-0-28 

— 

bz 

14-70 

— 

Zo 

0 

0 

Mo 

0 

0 

bo 

0 

0 

The  pitching  moment  derivatives  have  been  adjusted  to  the  correct  centre  of  gravity  position 
for  each  flight.  The  parameters  relating  specifically  to  the  non-linearities  of  the  pitching  moment 
I curve,  namely  bi  and  ba  are  absent  for  flight  21 1 since  in  that  case  the  moment 

curve  is  reasonably  linear  over  the  incidence  range  of  interest. 

I The  a priori  value  of  ZqImVe  has  been  set  to  zero  since,  for  a delta  wing  it  is  expected  to 

i be  small  and  is  difficult  to  estimate.  In  later  calculations  ZyjmVe  was  neglected  as  it  was  found 

to  have  very  little  effect  on  the  fit  error  or  on  the  other  estimated  parameters.  For  ai,  oa,  bi. 
ba  and  ba  the  a priori  values  were  calculated  by  assuming  G = I and  A = 0-3  in  the  formulae 
summarised  by  Equation  21 . The  parameters  ba  and  ba  were  only  used  in  later  runs  and  previously 
I identified  values  of  and  were  used  to  calculate  the  a priori  values  shown  in  Table  2. 

I The  elements  of  Da  fix  the  weighting  to  be  assigned  to  each  of  the  a priori  values.  These 

I were  set  by  assuming  a standard  deviation,  a,  equal  to  20%  of  the  parameter  value.  The  relevant 

element  of  Da  is  then  equal  to  \/a^.  Very  little  difference  in  the  identified  results  was  obtained 
by  assuming  a 10%  standard  deviation  but  smaller  values  of  a tend  to  push  the  identified 
parameter  values  towards  the  a priori  values,  especially  with  the  more  weakly  identified 
parameters.  This  effect,  together  with  the  relative  magnitudes  of  the  Cramer-Ra^  bound  helps 
^ to  isolate  the  weaker  parameters  (i.e.  those  that  only  weakly  influence  the  results)  if  they  are 

not  already  known.  Referring  to  Table  2,  those  parameters  with  zero  a priori  value  are  ascribed 
a large  standard  deviation  (99999  say)  which  leads  to  an  effectively  zero  weighting  element  in  Da. 

The  elements  of  Di  fix  the  weighting  to  be  assigned  to  each  of  the  time  histories  being 
matched,  e.g.  + a,,  q and  gotlVt  in  the  present  case.  No  attempt  was  made  to  relate  the 
elements  of  D\  to  the  instrument  accuracies  but  rather,  in  order  to  give  approximately  equal 
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weighting  to  each  matched  time  history,  a standard  deviation,  wm.  equal  to  5%  of  the  max  mum 
value  of  each  measured  record  was  selected  and  the  relevant  element  of  Di,  equal  to  l/w^,  then 
calculated.  The  resulting  weights  are  summarised  in  Table  3 below. 


TABLE  3 

Elements  of  Weighting  Matrix  Di 


Response 

Element 

Flight  210 

Flight  21 1 

a 

0 001 

0 001 

9 

5328 

6400 

0 

0 001 

0 001 

(am  + ««) 

91827 

81633 

gHilK 

147929 

49383 

For  records  not  being  matched  the  corresponding  element  of  Z)i  was  set  to  a small  value. 
This  applies  to  a and  0 in  Table  3.  In  the  case  of  0,  no  measurements  were  available  for  matching 
while  in  the  case  of  a,  the  corrected  value  + a,  is  being  matched  rather  than  a itself. 


4,  RESULTS 

In  this  section  the  results  of  the  identification  are  presented.  Two  particular  aspects  which 
arose  during  the  matching  process  are  first  discussed  in  some  detail.  One  involves  the  angle  of 
incidence  vane  corrections,  illustrated  by  reference  to  flight  211,  while  the  other  describes  the 
influence  of  non-linearity  in  the  pitching  moment  vj.  incidence  curve,  prominent  only  in  flight  210 
at  M = 0-96.  Following  this,  some  general  observations  arising  out  of  the  results  from  both 
flights  will  be  made. 

4.1  Matching  of  Angle  of  Incidence  Measurements 

Table  4 summarises  the  results  obtained  from  flight  21 1 and,  as  a reference,  shows  in  the 
first  column  the  a priori  values  and  weights  (based  on  20%  of  a priori  value)  of  the  various 
parameters  of  interest. 

The  column  headed  21  l-l  lists  the  results  obtained  when  only  the  qm  and  records  were 
matched  using  the  basic  mathematical  model  described  in  Section  2.1.  The  resulting  fits  to  qm 
and  were  very  good  as  shown  by  the  small  values  of  the  mean  square  weighted  errors  s\q) 
and  sHgtizIVe)  obtained.  The  identified  values  of  ZJmVe,  Z/JniVe  and  M)JI„  differ  only  slightly 
from  the  a priori  values,  while  MJfy  is  approximately  1 1 % down.  The  largest  change,  however 
is  in  the  damping  parameter  (A/j  + M^jly  which  is  little  more  than  half  the  a priori  value. 

Although  for  case  211-1  no  attempt  was  made  to  match  the  am  record  it  is,  nevertheless, 
of  interest  to  compare  the  mathematical  model  calculation  of  incidence,  a,  with  the  measured 
value,  am-  This  is  shown  in  Figure  4a  and  is  obviously  not  very  satisfactory.  If,  now,  the  uncor- 
rected am  record  is  included  in  the  matching  process,  together  with  qm  and  the  results  of 
column  211-2  apply  and  the  match  achieved,  for  angle  of  incidence  is  shown  in  Figure  4h. 
This  is  an  improvement  on  the  results  of  Figure  4a  but  is  far  from  being  good  with  a relatively 
large  mean  square  error,  5*(a)  of  2-78.  In  addition,  the  match  achieved  with  the  qm  and 
records  has  deteriorated  markedly  compared  to  case  21  l-l  with  sH,q)  going  from  0- 18  to  114 
and  .i^gritlVe)  going  from  0 078  to  0- 14.  The  values  of  the  identified  derivatives  also  vary  from 
that  obtained  in  case  21 1-1.  Most  significantly,  the  attempt  to  reconcile  the  three  traces  am,  qm 
and  rtt^  has  led  to  quite  large  changes  in  the  values  of  Mtjly  and  (Af,  ^ 
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At  this  stage  no  corrections  have  been  applied  to  measured  incidence,  The  next  case, 
21 1-3,  corrects  for  pitch  rate  effect,  a,  (see  Section  3.2)  and  in  addition  augments  the  basic 
mathematical  model  (Section  2.1)  with  the  lag  equation  as  described  in  Section  2.4.  This  resulted 
in  considerable  improvement  in  the  mean  square  errors,  r*,  for  all  the  records,  as  shown  in 
Table  4,  and  also  produced  the  much  improved  angle  of  incidence  match  of  Figure  4c.  The  lag 
time  constant,  t,  was  identified  as  being  very  small  (less  than  0 018)  and  in  fact  much  of  the 
apparent  lag  of  Figure  4b  was  removed  by  the  pitch  rate  correction,  The  improvement  in 
match  has  been  achieved  largely  as  a result  of  the  low  identified  gain  term  (C  = 0-80)  which 
would  correspond  to  a calibration  slope  decrease  of  20%.  As  one  might  hope,  the  main  identified 
aerodynamic  derivatives  have  moved  closer  to  the  results  of  case  211-1  which  neglected  the  a„ 
record  in  the  matching  process. 

On  closer  examination  of  Figure  4c  it  is  apparent  that  the  lag  model  fails  to  reproduce 
accurately  the  initial  rise  in  the  measured  record.  This  inadequacy  eventually  led  to  the  abandon- 
ment of  the  lag  model  and  its  replacement  by  a model  which  corrects  for  the  dynamic  imbalance 
of  the  incidence  vane  as  described  in  Section  2.3.  This  model  was  used  to  match  the  measured 
records  of  (corrected  for  pitch  rate  effect,  a,)  qm  and  nz^  and  the  results  are  shown  in  column 
211-4  of  Table  4 and  in  Figure  5.  The  improved  fit  for  each  of  the  three  records  is  evident. 
The  improvement  in  the  angle  of  incidence  match  is  the  most  marked  with  s\oi)  decreasing  from 
0 -393  for  case  21 1-3  to  0 - 174.  Figure  5 shows  an  excellent  match  for  all  the  response  variables 
and  in  particular  the  initial  portion  of  the  incidence  curve  is  now  accurately  matched.  In 
addition,  the  values  of  the  identified  aerodynamic  derivatives  are  now  very  similar  to  those 
obtained  from  case  211-1. 

Using  these  values  of  the  derivatives  and  the  identified  values  of  <22  and  bi  an  estimate  can 
be  made  for  K from  the  definitions  of  ai  and  b\  in  Equation  21.  This  results  in  values  of 
A = 0-27  and  0-31  from  a^,  bi  respectively  compared  with  the  value  of  A =0-30  used  in 
calculating  a priori  values  for  <22  and  bi.  However  neither  222  nor  bi  is  a strongly  determined 
parameter  as  is  apparent  from  the  small  changes  obtained  in  the  identified  values  and  Cramer- 
Rao  bounds  compared  with  the  a priori  values  and  weights.  On  the  other  hand,  ai  is  a more 
strongly  determined  parameter  and  using  a value  of  A = 0-3  together  with  the  estimated  value 
for  <21  implies  a value  for  G of  1 -01,  which  suggests  that  the  vane  calibration  is  indeed  close  to 
the  wind  tunnel  for  this  case. 

4.2  Non-linear  Pitching  Moment  Curve 

The  results  from  flight  210  are  summarised  in  Table  5 with  the  a priori  values  shown  in  the 
first  column  for  reference. 

Initially  only  the  qm  and  nz^  records  were  matched  using  the  basic  mathematical  model  of 
Section  2.1.  The  results  are  shown  in  the  column  under  210-1  in  Table  5 and  the  matched  q 
record  shown  in  Figure  6a.  The  match  is  nowhere  as  good  as  obtained  in  the  comparable  case 
211-1  of  the  previous  section.  On  closer  examination  of  Figure  6a  it  can  be  seen  that  while  in 
the  early  part  of  the  record  the  measured  value  of  qm  appears  to  have  a phase  lead  over  the 
calculated  q,  the  reverse  applied  in  the  latter  part  of  the  record.  This  systematic  phase  shift  is 
clearly  illustrated  in  Figure  66  which  plots  the  difference  between  measured  and  calculated  values, 
or  the  residuals,  against  time. 

The  match  can  be  greatly  improved  by  accounting  for  the  non-linearities  in  the  pitching 
moment  and/or  lift  curve  with  incidence  as  outlined  in  Section  2.2.  With  the  derivatives 
A/.i  and  Z.t  added  to  the  model  and  the  pitch  rate,  qm,  and  normal  acceleration,  nz^,  records 
again  matched,  the  results  are  as  shown  under  column  210-2  in  Table  S and  a comparison  of 
measured  and  calculated  pitch  rate  is  shown  in  Figure  7a.  In  calculating  the  and  elements 
of  the  input  vector  (Equation  22),  the  measured  angle.  corrected  for  pitch  rate  effect,  a«, 
was  used  in  this  case.  The  improvement  achieved  by  this  exter.d»l  model  is  quite  marked  with 
mean  square  error  s\q)  decreasing  from  0-459  to  0- 169  aiid  sHgn,iy,)  decreasing  from  0-393 
to  0-169.  The  plot  of  the  residuals  in  Figure  76  also  illustrates  the  improved  match.  For  an  ideal 
match  the  residual  at  any  one  point  should  be  uncorrclated  with  that  of  any  other  point  on  the 
record.  To  check  this  the  autocovariance  function  for  the  ^residuals  was  calculated.  This  is 
defined  as  j. 

Aq)  = + t)/^  (30) 
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where  the  bars  refer  to  the  time  average  value  and  t is  a given  time  shift.  The  plots  ot f{q) 
against  t for  the  linear  and  non-linear  models  are  shown  in  Figure  8 and  confirm  the  improvement 
achieved  using  the  non-linear  model.  Although  the  pitch  rate  records  have  been  used  to  illustrate 
the  points  made  so  far,  the  same  conclusions  apply  to  the  normal  acceleration,  rti,  records. 
It  may  be  noted,  finally,  in  comparing  cases  210-1  and  210-2  that  although  the  match  has  been 
significantly  improved,  the  estimated  values  of  the  aerodynamic  derivatives  have  altered  only 
slightly  with  the  exception,  as  expected,  of  MJIy. 

The  next  case,  210-3,  matches  all  three  records,  a™,  qm,  and  The  non-linear  model  for 
the  pitching  moment  curve  was  retained  (i.e.  the  M^i,  A/,j  terms)  but  a linear  model  for  the 
Z-force  curve  was  assumed  (i.e.  the  2,2  term  was  removed)  since  the  improvement  achieved  by 
retaining  the  term  was  relatively  small.  The  angle  of  incidence,  a^,  was  corrected  for  pitch 
rate  effect,  aj,  and  the  corrected  value,  otw  -f  ««  used  in  calculating  the  a*  and  a®  terms  of  the 
extended  input  vector.  The  complete  mathematical  model  also  included  corrections  for  the 
dynamic  imbalance  of  the  incidence  vane  as  described  in  Section  2.3,  which  led  to  the  additional 
parameters  ai,  aa,  b\,  bi  and  63.  The  results  of  case  210-3  are  summarised  in  Table  5 and  plots 
of  all  relevant  variables  shown  in  Figure  9.  An  excellent  match  is  obtained  in  all  cases.  The 
relatively  large  mean  square  error  in  the  incidence  match,  i^fa),  reflects  the  higher  noise  levels 
on  the  incidence  measurements  compared  to  the  other  measured  quantities  (see  Figure  9).  The 
slight  increase  in  the  value  of  sHgUzIVe)  compared  to  case  210-2  can  be  ascribed  to  the  neglect 
of  Za2.  The  values  of  the  identified  parameters  are  very  similar  to  those  from  case  210-2.  The 
values  of  K implied  by  at  and  as  are  0-30  and  0-32  respectively  while  the  estimated  value  of  ai 
implies  a value  of  C = 0-84  (for  K = 0-30).  This  low  value  of  G suggests  a substantial  change 
in  the  incidence  vane  calibration  curve  slope  for  this  case  (M  = 0-96)  in  contrast  to  results  in 
the  previous  section  from  flight  21 1 (M  = 0-71). 

The  final  case,  210-4,  is  a repeat  of  210-3  with  the  sole  difference  that  instead  of  using  the 
measured  values  of  incidence  («m  + «»)  to  calculate  the  1*2  and  a*  terms  of  the  input  vector, 
these  terms  were  calculated  using  the  values  of  the  state  variable  a as  obtained  in  210-3.  The 
results  in  Table  5 show  only  minor  changes  from  the  results  of  case  210-3  with  the  exception 
of  the  extracted  values  of  A/5,2  and  A/»j. 

The  significance  of  the  values  of  A/5,2  and  can  be  illustrated  by  reference  to  Figure  10, 
which  plots  the  pitching  moment  against  total  incidence,  ««  f a.  The  identified  pitching  moment 
takes  the  form 

A/  = A/(a  = Xf)  -f-  AM 

where  the  incremental  moment,  AM,  is  (according  to  Equation  13); 

AM  = M^x  + M^rx-’-  4-  A/5„a3 

and  M{x  = a«)  is  taken  from  the  wind  tunnel  results*  of  pitching  moment  v.y.  incidence,  adjusted 
to  the  correct  centre  of  gravity  position.  The  values  of  A/5,2  and  M^y  from  cases  210-2  and  210-3 
produce  almost  identical  curves  over  the  incidence  range  covered  in  the  pitch  up  manoeuvre 
despite  the  differences  in  the  actual  values  of  A/„2  and  M^i.  It  is  clear  that  a fairly  wide  range 
of  values  of  A/5,2  and  A/^j  can  lead  to  almost  identical  curves  over  the  restricted  range  of  interest 
and  this  accounts  for  the  large  Cramer-Rao  bounds  obtained  for  these  parameters.  The  curve 
corresponding  to  case  210-4  is  also  shown  in  Figure  10  and  differs  slightly  from  the  other 
calculated  curve  towards  the  extremes  of  the  a-range.  For  comparison,  the  equivalent  wind 
tunnel  data®  at  A/  = 0-95,  is  also  plotted  on  Figure  10  and  appears  to  agree  very  well  with  the 
identified  results.  This  is  strong  support  for  the  approach  used  and  confirms  that  non-linear  aero- 
dynamics can  be  accommodated  successfully  within  the  framework  of  a linear  mathematical  model. 

4.3  General  Discussion 

In  the  present  calculations  the  Cramer-Rao  bounds  are  not  strictly  correct  since  they  are 
based  on  the  prescribed  D\  and  Dt  matrices  rather  than  the  covariance  matrices  as  required 
by  maximum  likelihood  theory  (see  References  5 and  6).  Nevertheless  the  relative  magnitudes  of 
the  calculated  bounds  are  an  indication  of  the  amount  of  confidence  to  be  placed  in  the  identified 
results.  This  is  particularly  true  when  the  a priori  feature  is  used.  In  that  case,  if  the  bounds  do 
not  differ  from  the  a priori  weights  prescribed,  it  is  an  indication  that  little  additional  information 
is  contained  in  the  data  and  the  a priori  value  should  be  accepted  for  that  particular  parameter. 
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Hence,  looking  at  the  final  results  of  Flight  21 1 (case  211-4)  and  210  (cases  210-3  and  210-4) 
in  Tables  4 and  5,  it  is  apparent  that  is  a weakly  identified  derivative  and  the  a priori  value 
should  not  be  modified  on  the  basis  of  the  test  information  available.  On  the  other  hand, 
reasonable  confidence  can  be  placed  in  the  identified  values  of  the  other  four  derivatives,  ZJmVe, 

For  flight  211,  ZJmVe  is  in  good  agreement  with  the  a priori  value  but  for  flight  210,  at 
M = 0-96,  it  is  significantly  less,  being  0-74  compared  to  0-85.  The  rapid  flow  changes  which 
occur  at  transonic  speeds  may  be  the  cause  of  this  difference.  The  estimated  values  of  MJIy 
and  M)Jly  from  flight  211  are  both  approximately  10%  less  in  magnitude  than  the  a priori 
values  (see  footnote  to  Table  4).  This  sort  of  difference  could  easily  be  accounted  for  by  an 
error  in  centre  of  gravity  position  since  the  corrections  for  shift  from  the  wind  tunnel  reference 
position  (52%)  are  quite  large  in  this  case.  Alternatively,  a 10%  error  in  the  pitch  moment  of 
inertia,  ly,  could  also  account  for  the  difference. 

For  flight  2l0  good  agreement  is  obtained  with  Mgjly  but  MJIy  appears  to  be  considerably 
lower  than  the  a priori  value.  However,  as  discussed  in  Section  4.2,  the  slope  of  the  M v.s.  a curve 
varies  considerably  at  A/  = 0-96  and  inspection  of  Figure  10  shows  that  the  slop)e  can  alter  by  , 
a factor  of  2 over  the  range  of  interest.  In  fact.  Figure  10  indicates  a reasonably  good  agreement 
between  the  estimated  and  wind  tunnel  pitching  moment  curves. 

Looking  now  at  the  damping  derivative,  (Mq  + Mi)lly,  the  estimated  value  at  A/  = 0-71 
(flight  21 1)  is  little  more  than  half  the  a priori  value  while  at  A/  = 0-96  (flight  210)  it  is  in  good 
agreement  with  the  a priori  value.  The  low  value  at  M = 0-71  cannot  be  easily  explained  away 
in  view  of  the  relatively  good  agreement  of  the  other  moment  derivatives  and  the  conclusion 
is  that  the  a priori  value  is  too  high  and  the  lower  value  should  be  used  in  any  future  mathematical 
model.  The  Datcom  estimate* ^ gives  a value  of  0-94  for  (A/,  t-  M^jly  at  Af  = 0-7  with  the 
c.g.  at  the  49- 1"'„  position.  This  is  almost  25%  below  the  present  estimated  value.  Further  the 
value  for  Mqjly  of  about  0-68  quoted  for  the  Mirage  III  in  Reference  12  is  also  in  line  with  the 
Datcom  estimates.  These  values  confirm  the  conclusion  that  the  a priori  value  used  in  the  present 
case  at  A/  = 0-71  is  too  high.  On  the  other  hand  the  good  agreement  at  A/  = 0-96  suggests  that 
the  shape  of  the  damping  curve  as  Mach  number  changes  is  different  to  what  had  previously 
been  assumed.  Further  work  is  planned  to  establish  the  variation  of  damping  over  the  range  of 
Mach  numbers  from  0-7  to  1-2. 

Some  comment  should  be  made  on  the  choice  of  records  for  the  matching  process.  It  has 
been  possible  to  achieve  excellent  matches  for  all  the  available  records,  i.e.  q„  and 
provided  the  incidence  vane  errors  were  properly  accounted  for.  It  has  also  been  shown  that 
successful  identification  can  be  achieved  using  qm  and  records  only.  For  example,  com- 
parison of  211-1  with  211-4  or  210-2  with  210-4  show  close  agreement  in  the  values  of  the 
estimated  aerodynamic  derivatives  and  equally  good  matches  of  the  respective  time  response 
records.  Thus  if  doubt  exists  about  the  accuracy  of  a particular  record,  exclusion  of  that  record 
can  still  lead  to  good  results.  If  however,  an  attempt  is  made  to  match  an  erroneous  record 
e.g.  uncorrected  incidence  measurements,  then  an  all  round  deterioration  in  results  can  occur. 
For  example,  compare  case  2II-I  with  211-2.  Similarly,  an  inferior  result  is  obtained  if  non- 
linearities  in  the  M vs.  « curve  are  not  modelled  e.g.  compare  case  710-1  with  210-2. 


5.  CONCLUDING  REMARKS 

A modified  Newton- Raphson  parameter  estimation  procedure  has  been  applied  to  flight 
data  from  a 60°  delta  wing  aircraft.  The  response  of  the  aircraft  over  a period  of  5 seconds 
(ou„r=  21)  to  an  elevon  pulse  input  at  two  Mach  numbers,  0-71  and  0-96,  was  successfully 
analysed  to  obtain  all  the  relevant  longitudinal  aerodynamic  parameters.  The  sampling  rate  of 
60  per  second  (wnAr  — 0-07)  was  found  to  be  ample  and  no  problems  of  numerical  convergence 
were  encountered.  Records  of  elevator  input  and  the  response  variables  of  incidence,  pitch  rate 
and  normal  acceleration  were  used  in  the  analysis.  Good  results  were  achieved  when  only  pitch 
rate  and  normal  acceleration  were  matched  and  this  was  confirm--^  by  matching  all  three 
variables  simultaneously. 

In  order  to  match  the  incidence  measurements  satisfactorily  it  v'>s  found  necessary  to 
correct  for  pitch  rate  induced  error  on  the  incidence  vane  and,  in  addition,  to  make  allowance 
for  inertial  contributions  to  the  vane  readings,  which  were  not  dynamically  balanced,  and  for 
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possible  wind  tunnel  calibration  errors.  The  vane  calibration  factor  which  produces  the  best 
match  is  one  of  the  parameters  obtained  from  the  program  output  and  this  can  then  be  checked 
against  the  wind  tunnel  calibration. 

A technique  for  including  non-linear  aerodynamic  characteristics  within  the  context  of  a 
linear  model  has  been  demonstrated.  The  approach  used  produced  a pitching  moment  curve  in 
good  agreement  with  wind  tunnel  test  results.  Making  allowance  for  the  non-linearity  of  the 
pitching  moment  curve  also  achieved  significant  improvement  in  the  matching  of  the  pitch  rate 
record  at  A/  = 0-96.  The  present  method  of  analysis  could,  with  possible  advantage,  also  be 
applied  to  extraction  of  data  from  wind  tunnel  tests  with  oscillating  models. 

The  estimated  values  of  the  longitudinal  derivatives  obtained  in  the  present  study  seem  to 
be  in  reasonable  agreement  with  wind  tunnel  values.  One  exception  appears  to  be  damping  in 
pitch  where  the  results  disagree  with  the  a priori  value  at  M = 0-71  but  are  in  good  agreement 
at  Af  = 0-96.  This  important  area  is  being  studied  further  with  the  aim  of  establishing  the 
variation  of  the  pitch  damping  over  the  Mach  number  range  from  0-7  to  1-2. 

An  important  aim  of  the  work  reported  here  has  been  to  test  the  application  of  systems 
identification  procedures  to  available  longitudinal  flight  data.  The  results  have  demonstrated  a 
variety  of  ways  in  which  the  program  can  be  used,  not  only  to  validate  mathematical  models 
and  check  flight  results  against  data  obtained  from  wind  tunnels,  but  also  to  obtain  information 
not  easily  available  through  other  means. 
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Coefllcienls  in  Equation  21 

Matrix  of  stability  parameters.  Equation  2 

Coefficients  in  Equation  21 

Measurement  bias  vector 

Equation  error  bias  term.  Equation  23 

Matrix  of  control  parameters.  Equation  2 

Unknown  parameter  vector 

Vector  of  a priori  parameter  values 

Mean  aerodynamic  chord 

Non-dimensional  moment  coefficient,  Mj  1/2  p V^Sc 
Weighting  matrix  for  response  variables.  Equation  I 
Weighting  matrix  for  parameter  estim:  es.  Equation  I 
Matrix  of  stability  parameters.  Equation  3 
Gravitational  acceleration 

Matrix  of  control  parameters.  Equation  3 or  gain  of  vane.  Equation  17 
Height  above  sea  level,  m 
Moment  of  inertia  in  pitch,  kg-m* 

Cost  functional  (Fit  error) 

Constant  of  proportionality.  Equation  18 
mass,  kg 

Mach  number  or  pitching  moment 
Equation  error  bias  term.  Equation  10 
Pitching  moment  derivative  w.r.t.  q 
Pitching  moment  derivative  w.r.t.  u 
Pitching  moment  derivative  w.r.t.  a 
Pitching  moment  derivative  w.r.t. 

Pitching  moment  derivative  w.r.t.  a® 

Pitching  moment  derivative  w.r.t.  i 
Pitching  moment  derivative  w.r.t.  S* 

Measurement  noise  vector 
Normal  acceleration  in  g units 
Number  of  data  samples 
Roll  rate,  rad/s 
Pitch  rate,  rad/s 
Yaw  rate,  rad/s 
Mean  square  weighted  error 
Reference  wing  area,  m^ 

Time,  sec 

Length  of  measured  responses,  sec 

Vector  of  control  input  variables 

Velocity  increment  in  x-direction,  m/s 

Resultant  airspeed,  m/s 

Velocity  increment  in  2-direction,  m/s 

Vector  of  state  variables 

Body  axis  co-ordinate  in  forward  direction 

Force  in  x-direction,  N 

Y-force  derivative  w.r.t.  u 

Y-force  derivative  w.r.t.  a 


X-force  derivative  w.r.t.  Se 

y 

Calculated  response  vector,  Equation  3 

y 

Body  axis  co-ordinate  in  lateral  direction 

z 

Measured  response  vector,  Equation  4 

7 

Body  axis  co-ordinate  in  down  direction 

z 

Force  in  7-direc(ion,  JV 

Zo 

Equation  error  bias  term.  Equation  10 

Zu 

Z-force  derivative  w.r.t.  u 

z. 

Z-force  derivative  w.r.t.  q 

Z. 

Z-force  derivative  w.r.t.  « 

Z.z 

Z-force  derivative  w.r.t.  ot^ 

2.3 

Z-force  derivative  w.r.t. 

Z-force  derivative  w.r.t.  8, 

a 

Angle  of  incidence  increment,  rad 

ajT 

Pitch  rate  correction.  Equation  16 

Inertial  correction  to  «,  Equation  18 

A 

Increment 

S. 

0 

Elevator  angle  increment,  positive  t.e.  down,  rad 
Pitch  attitude  increment,  rad 

P 

Air  density,  kg/m® 

a 

Standard  deviation 

<^CR 

Cramer-Rao  bound 

T 

Time  constant.  Equation  25.  or  Time  shift.  Equation 

W„ 

Natural  undamped  frequency,  rad/s 

Subscripts 

a 

Normal  accelerometer 

e 

Trim  or  Equilibrium  state 

i 

Time  index 

m 

Measured  value 

V 

Incidence  vane 
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FIG.  1 BASIC  CONCEPT  OF  ESTIMATION  METHOD 


FIG.  2 SMOOTHED  ELEVATOR  INPUT  - FLIGHT  21 1 
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FIG  3 SMOOTHED  ELEVATOR  INPUT  - FLIGHT  210 


\d.ez 


t — seconds 

(a)  Basic  model  and  measured  incidence  (unmatched) 


FIG.  4 COMPARISONOF  ANGLE  OF  ATTACK  FOR  FLIGHT  21 1,  M = 0.71 
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FIG.  5 FINAL  MATCHED  RESULTS  FOR  FLIGHT  21 1,  M = 0.71 
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FIG.  6 COMPARISON  OF  PITCH  RATE  FOR  FLIGHT  210,  M ■=  0.96,  LINEAR  MODEL 
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FIG.  7 COMPARISON  OF  PITCH  RATE  FOR  FLIGHT  210,  M = 0.96,  NON-LINEAR  MODEL 
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FIG.  9 FINAL  MATCHED  RESULTS  FOR  FLIGHT  210,  M 0.96 
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FIG.  10  COMPARISON  OF  WIND  TUNNEL  AND  IDENTIFIED  Vs  a CURVES 
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